b 



Two-dimension plasma expansions with anisotropic pressure 



Yongsheng Huang/- ^ G Yuanjie Bi/-^ Xiaojiao Duan/ 
Naiyan Wang/ Xiuzhang Tang/ and Zhe Gao^ 

'China Institute of Atomic Energy, Beijing 102413, China. 
^Department of Engineering Physics, Tsinghua University, Beijing 100084, China. 

(Dated: November 28, 2008) 

Abstract 

A two-dimension self-similar solution is proposed for a plasma expansion with anisotropic pressure. 
With the solution, it depends on the relationship between the ratio of the longitudinal and the transverse 
temperature of the plasma, and the electron-ion mass ratio, jj., that the plasma front is composed by a part 
of hyperbolic (or a plane) and a small pointed projection at the center or a part of an ellipse. Zhang and 
coworkers 's experiments (PRL, 99, 167602 (2007))support our results for /c^ € (t, 1]. For <t, there is an 
anomalous high-energy plasma emission at the angle of near 90° due to longitudinal Coulomb explosion. 
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FIG. 1: (Color online) (^1,^2) at the ion front for different k^, the ratio of Tj_ and Ty, which is compared 
with n = Zmelnii, the mass ratio of electron and ion for cjp - \ and Co = in Eq. Q. Figure [Ha) and 
(b) are obtained from the solutions of Eq. ([2]). In Figure (Ha), = Tj_/T\\ > fx. In these cases, the curve 
of the ion front is a part of an ellipse, especially, a cycle for = I. The major axis of the ellipse is in the 
lower-temperature direction. In Figure [T](b), = T±/T\\ < fi. In this case, the ion front contains two parts. 
For ^2 ^ 2.5, the ion front is a part of a hyperbola or a plane (k^ = fj.). For ^2 ^ 2.5, the ion front is a small 
pointed projection. 

The picosecond and femtosecond laser ablation yj], the ultrafast diagnostics of hydrodynam- 
ics M\, the laser-ion acceleration [3], etc., attract more and more international attention since the 
great development of ultrashort and ultraintense laser technologies. The theory of plasma ex- 
pansions is not only the fundamental mechanism of laser ablation but also the laser-ion acceler- 
ation. Although various theory such as phase explosion [4], Coulomb explosion [5^, self-similar 
plasma expansions 1^, [vl], two-dimension self- similar plasma expansions |8], and a relativistic 
model[^ have been proposed, the observed angular-ion distributions such as ring structure [fioll . 
time-resolved elliptic shadowgraphs of material ejection [1] and coronal hydrodynamics of laser- 
produced plasma[2] have not been explained analytically. The problems are that: the electrons 
generated in the laser-solid interactions are anisotropic generally, however the self-similar expan- 
sion JgI Q] and the relativistic model are both one-dimension, the two-dimension theory is only 
valid for a plasma with isotropic pressure, therefore, they can not predict the shapes of plasma 
front with anisotropic pressure. 
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In this letter, a self similar two-dimension solution for a plasma expansion with anisotropic 
pressure is obtained. Although isothermal-expansion assumption is used, the solution is still suit- 
able for a system whose temperature changes slowly, since the whole process can be separated to 
several periods and the temperature can be assumed as a constant in any period. In the solution, 
there are two important parameters: the ratio of the longitudinal temperature, and the transverse 
temperature, T\\ of the plasma, = Tj_/T\\, and the electron-ion mass ratio, fi = Znie/nii, where Z 
is the charge number of the ion, nieinii) is the electron (ion) mass. With the solution, it is found 
that the relationship between and decides the shape of the ion front: a complicated surface or 
an ellipse shown by Figure [IJ For 6 [0,/z), the ion front is composed by a part of a hyperbolic 
and a small pointed projection at the center. In the critical case, = fi, the ion front is a plane 
and a small pointed projection at the center. For e the ion front is a part of an ellipse and 

the major axis is in the lower-temperature axis. The ion velocity at the ion front has the similar 
construction to that of the ion front. However, for 6 (0, 1), the major axis of the ion- velocity 
ellipse is in the higher-temperature axis. The difference of the angular-energy distribution from the 
known onefsi H ] is that the energy is a delta function at the maximum angle of near 90° for < jj.. 
It is an anomalous high-energy plasma emission at the angle of near 90° due to the longitudinal 
Coulomb explosion. The plasma expansion is adiabatic and the temperature changes slowly in a 
vacuum. However, there is energy transmission from the plasma to the air besides that between the 
two directions of the plasma. With the assumption: the plasma temperature decreases with time 



exponentially in the air[|12l]. the dependence of the position of the ion front in the central axis on 
time is compared with that observed by Zhang and coworkers 's in [1]. From the comparison, it is 
inferred that the plasma front shown in [jl|] should be the front of Al^, Al^^ or AP^. 



We assume the pressure tensor of a two-dimension anisotropic plasma satisfies : 

(1) 



neksT^ 

if the ion temperature hgTi «; ksTw and hgTi ksT^, where T\\ is the electron temperature in the 
x direction and is the electron temperature in the y direction which is perpendicular to x. In the 
following discussion, we will show the importance of the ratio, Tj_/T\\, and the dependence of the 
shape of the plasma front on the ratio. In the calculation of the formula of the ion front, we assume 
the plasma expansion is isothermal. However, if the temperature changes with time slowly, in any 
interval short enough, the temperature can also be considered as a constant and our calculation and 
results still hold. Therefore, to a certain extent, our calculation and results allow a wide latitude of 



generality in applications. 

For convenience, the physical parameters: the time, t, the length coordinate, x{y), the ion (elec- 
tron) velocity, V;(Ve), the electron field, E, and the ion (electron) density, n,(ne) are normalized by 



the inverse of the equivalent plasma frequency, cOpiQ = ^fn^^e^^Jm^, the equivalent plasma De- 
bye length, AuQ = cjcopio, the equivalent ion acoustic speed, = yfZk^TjTni, Eq = kBTe/eAoo, 
and the reference hot-electron density, n^o respectively, where m, is the ion mass, Z is the charge 
number of the ion, e is the elemental charge and Te = T\\ + is an equivalent temperature. 
Then the electric potential is normalized as = eij/lksTe, where ij/ is the physical potential. 
t,x{y),Vi{Ve),E,ni{ne) are still used to represent the normalized parameters in the following dis- 
cussion. 

In the new frame: t = t,^\ = x/Ri{t),^2 = y/Riit) , the ion (electron) velocity[8] sat- 
isfies: Vi^x(Ve,x) = ^i^p v,;y(Vg,j,) = ^2^2- "'"he iou (electron) density is assumed as: = 
Ni{e),i(^i,^2)/Rl+Ni(e),2(^i,^2)/Rl, wMch is similar to that in [|8D. With the transformation, the equa- 
tions of continuity and motion, and Poisson's equation are obtained easily in the new coordinate 
system. The condition for the automatical satisfiability of the continuity equation is R2/R1 = k, 
where a- is a constant. Then the ion (electron) density is simplified to Ni(e){^\,^2)IR^{t), where 
Ni(e) = + Ni(ei2/K^, R = Ri- Solviug the new motion equation of ions gives the potential in 
the ion region, (p = -(p^i^j + k^^I), where 0° is a constant. Therefore the electron density in the ion 
region satisfies: Ne{^\,^2) = A^e,o exp[-(l + ix)(p^{^\ + ^\)la\\ by solving the electron motion equa- 
tion with = Tj_/T\\, where ai = T\\/Te. From Poisson's equation, the ion density is A^, = Ne + 4. 
The electron motion equation can be solved in the ion region due to the special form of the electric 
potential: (p = 0i(^i) + 02(^2)- However, beyond the ion front, the potential and electron density 
are governed by the motion equation of electrons and Poisson's equation together, since the po- 
tential can not be separated to 0i(f 1) -1- ^2(^2) and then the electron density can not be solved from 
the motion equation solely. Combined the motion equation and Poisson's equation, a two-order 
partial difl'erential equation of the electron density can be achieved. The first integral of it gives: 
+ /c2(g)2 = 2^'^p(y)"4(i+^^Wt- ^ Q -g ^j^g integral constant and Y = \n(N,). 

The physical condition requires that Y is a function and therefore, the curve equation of the ion 
front is: 

A 2 = 1 - HK^,B ^ = - H,D = 2(i+^)(iV;,2)^o.2 ■ 

where Co is a constant and cfP'^ is the square of 0°. With Eq. (O, the curve of the ion front may be 
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some part of a hyperbola, an ellipse or a cycle approximately for different ratio of the temperatures 
( Tx and ry), k^. A critical value of is the mass ratio of electron and ion, // or the inverse of it. 
Due to the assumed symmetrical form of the pressure, we just need to consider the cases for 
= Tj_/T\\ e [0, 1]. Figure [U shows the eight cases of the ion front for e [0, 1], which are 
obtained by solving Eq. @ with Cq = and 0" = 1. For > ju, the curve of the ion front is a 
part of an ellipse whose eccentricity is e = ^^ iJ^i^'^^ as shown by Fig. [Ua). As ^ ju, e ^ 1. 
Specially, for /c^ = 1, i.e., = Ty, the ion front is a part of a cycle. For < /i, the ion front 
is complicated and contains two parts: a part of a hyperbola and another curve as shown by Fig. 
[IJb). For the part of the hyperbola, ^2 > 2.5 and the slope of the asymptote is -yj^^^- As fx, 
the slope trends to infinite and the ion front becomes a flat surface as shown by the dash dot line 
inFig.mb). 

Now imaging the following physical picture: the initial state is T^.o = juT'ilo, and the temperature 
Tx increases slowly to the equilibrium temperature, Teq ( Note: Teq = (ju + l)r||,o/2 in our two 
dimensional case. However, in the three dimensional case, if we assume the temperature in the y 
direction is equal to that in the z direction at the very beginning time, T^q = (2// + l)r||_o/3. ), when 
the thermodynamic equilibrium is reached, since the energy is transmitted between the transverse 
direction x and the longitudinal directions (Note: in two dimensional case, y or z is the longitudinal 
direction; in the three dimensional case, y and z are all longitudinal directions where we assume 
the temperature in the y direction is equal to that in the z direction.) with collisions. At the very 
beginning time, the ion front is a flat surface with a small pointed projection at the center shown 
by the dash dot line in Fig. [IJb). After an interval, > /i, and then the ion front develops to a part 
of ellipses shown by the solid, dash, and dot line step by step in Fig. [2a), which are also shown 
by Figure 2 (b-h) in yj]. At last, the ion front trends to a part of a cycle shown by the dash dot line 
in Fig. [2 a), which is also shown by outer layer of Figure 2(i) in yJ]. That is a depiction of the two 
dimensional plasma expansion into a vacuum. However, it has not been verified experimentally in 
a vacuum. In spite of that, the time-resolved shadowgraphs of material ejection shown by Zhang 
and coworkers in dil] support our theory without doubt. The outermost layer of the shadowgraph 
does undergo the physical processes depicted by our above discussion. It is the most difference 
between the experiment and our scene that the plasma expansion is in the air or in a vacuum. In 
the air, there are thermal transmission from the plasma to the air, collisions between the plasma 
and the air molecules and other nonlinear effects, and then the outermost layer of the plasma cools 
down. In this case, the shock wave front (shown by Fig. 3 in yJ]) forms since the particles inner the 
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FIG. 2: (Color online) The ion velocity distribution at the ion front for / - lO^^W/cm^, A - 0.8nm, and 
^ Tj_/T\\ € [0, 1], 0° ^ 1 in Eq. ©, a{t) ^ 0.003 at t ^ O.lps. (a) The ion velocity distribution at the ion 
front for < jj.. (b) The ion velocity distribution at the ion front for > ji. 

plasma catch up with the particles at the front. In a vacuum, these phenomenons can not happen. 

With the ion motion equation, R{t) satisfies: Qx^{a\)da\ = ->/^t/RQ, where a = ^In^R/Ro) 
and Rq = R{t = 0). Figure |2] shows the ion-velocity distribution at the ion front obtained from 
our two-dimension theory for different temperature ratio, = T^IT\\ e [0, 1]. For different times, 
the shapes of the ion-velocity distribution are the same although the amplitudes are different. For 

6 (0,ju), the velocity distribution has two parts: ellipse-like at small V2 and hyperbolic-like at 
large V2, which is shown by the dot line or the dash line in Fig. Eta). The solid line in Fig. [2ta) 
shows all the ions move towards the x direction since the temperature in the y direction is zero. 
The dash dot line in Fig. [2ta) shows the critical case: = ju, in which the velocity in the x 
direction at large V2 is zero and all the ions move in the y direction. For /c^ 6 (//, 1], the ion velocity 
distribution is a part of an ellipse whose eccentricity is ^JuiKT^ - k^)I{1 - piK'^) and the major axis 
is in the high temperature direction, x. At /c^ = 1, the ion- velocity distribution is a part of a cycle 
and isotropic as shown by the dash dot line in Fig. [2tb). 

Figure |3] shows the angular-energy distribution for 6 [0, 1]. With it, the facts are indicated 
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FIG. 3: (Color online) The ion energy distribution at the ion front VS the angle to the central axis, x, for 
- Tj_/T\\ € [0, 1], 0*' = 1 in Eq. Q. In this figure, the energy for any angle is normalized by Eq for each 
K^, where Eq is the ion energy at the central axis. 

: (1) for e [0,//], the energy decreases with angle except for the maximum angle of near 90°. 
At the maximum angle, the energy is a delta function with the peak value of about 0.62£'o- Espe- 
cially, the maximum angle is 90" for the critical point, = /i. This anomalous plasma emission 
happens since the large Coulomb space-charge field dominates in the longitudinal direction and is 
counteracted partly by the gradient of the pressure in the transverse direction. (2) for 6 (/i, 1), 
the ion energy decreases with angle monotonously. At = 1, the energy is a constant function 
with respect to angle. 

As discussed above, there is energy transmission from the plasma to the air for plasma expan- 
sions in the air. Therefore, besides the energy transmission between the longitudinal direction and 
the transverse direction of the plasma, the temperature of the plasma decreases exponentially to 
zero for the time long enough since the temperature of the air is about OeV. As an important ap- 
plication, the position of the ion front in the air can b e g^ iven by our solution with the assumption: 
the temperature, decreases exponentially to zero[12] with the scale time of the temperature, 
teff = -d\n{Te)/dt. Therefore Te = Teo^^pi-t/teff) in the air and = (which is (r^o + T'no) 
in the two dimensional case and {2Tj_o + T\\o) in the three dimensional case where the longitudinal 
temperatures are the same.) in a vacuum. The ions at the outermost layer are decelerated for the 
plasma expansion in the air after a period of time shown by the dash line and the dash dot line in 
Figure Ufa), however, the velocity of them tends to a constant in a vacuum shown by the dash line 
and the dot line in Figure Ufa). 

Figure Sfa) shows different ion fronts in a vacuum, in the air and for different initial charge- 
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FIG. 4: (Color online) (a) The ion front of Al^+ in the central axis,x, with ^ij = 0.42,^2,/ = 0, VS t, for 
rieo = and Z?i,o = 26?ieo> shown by the solid line (in a vacuum) and the dash line (in the air, with 

the scale time of temperature, te/f = 2ns), and Z«;o = (1 + 2.5 x 10~')?ieo> shown by the dot line (in a 
vacuum) and the dash dot line (in the air, with the scale time of temperature, tgff - 2ns); (b) the front of 
different ions in the central axis ,x, with = 0.42,^2,/ - 0, VS t. The solid line, the dash line and the dot 
line correspond to the expansion in the air and are obtained by our solution. The solid line is for Al^, and 
Znio = 7.25?igo> teff - l-5ns. The dash hne is for AF^, and Z?i,o - 1.0625neo> teff - 2ns. The dot line is for 
Al^"^, and Zn,o = 2neo, teff - 2ns. The dash dot hne , the dash dot dot line, the short dash line correspond 
to Figure 4 in Uj]. 



separation densities. The difference of the ion fronts in a vacuum and in the air has been discussed 
in the above paragraph. In the air, the collisions between the plasma and the air molecules are 
important for the expansion and considered gracefully through the exponential decrease of the 
plasma temperature. The collisions result in nonlinear phenomenons and energy transmission. For 
different initial charge-separation densities of the plasma, the stronger the initial charge-separation 
field and the more efficient the ions are accelerated as shown by the comparison between the solid 
line and the dot line or between the dash line and the dash dot line in Fig. I3a). When an ultra- 
short and ultra-intense laser pulse interacts with a solid target, many interesting phenomenons: 
hot-electrons generation and acceleration, x-ray and y-ray emission, fast ions generation, plasma 
ablation and plasma splash and so on occur. For different phenomenon, the scale time is different. 
A fact is: the smaller the scale time, the stronger the charge-separation field and the nonlinearity 
in parameters of the plasma. For the generation of fast ions, the scale time is about 1 - lOps. 
However, for the laser-plasma ablation, the scale time is several nanoseconds or even longer. In the 
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period of several pulse durations, the charge-separation field is strong, i.e., the density difference, 
Zfii - He (which is here.) is large (since R{t) increases with time and tends to infinite as 
t — > +00.). Therefore, with the self-similar solution for energetic ions generation given by Huang 
and coworkers' in [|7D, the high charge-mass ratio, Z/M, ions are accelerated preferentially and 
efficiently. When an ultra-intense laser pulse interacts with a solid target, H^, C"^^, C^^, O"^^ and so 
on are accelerated orderly. For Al target, Al^^, Al^^, Al^ are accelerated lastly. According to our 
analysis, the front surface observed by Zhang and coworkers' in should be the ions of Al and 
the density of the ions should be low since the scale time is of nanoseconds. The fronts of the ions: 
Al^+, Al^+, A1+ in the central axis, X, with ^i j = 0.42,^2,/ = are calculated and compared with 



the curves given by Zhang and coworkers' in [10 as shown by Figure |4](b). They are consistent. 

In conclusion, the analytic solution for the ion front of plasma expansion with anisotropic pres- 
sure is predicted. It is supported by experiments |lli|20 lately. The ion- velocity distribution at the 
ion front, angular-energy distribution at the ion front are both given. It needs to be concerned that 
the energy is a delta function at the maximum angle of near 90° and the plasma emits anomalously 
due to the longitudinal coulomb explosion for e [0, //] as shown by the solid, dash, dot lines in 
Fig. [31 In a vacuum, the expansion is adiabatic, however, in the air, the temperature decreases with 
time exponentially and collisions happen, and then the ions at the outermost layer are decelerated 
and shock waves generate in the front of the plasma front. From the comparison (shown by Fig. 
131) , it is concluded that the ions observed in yj] are the ions of Al. However, the collisions and 
the nonsymmetric components of the pressure tensor of plasmas have not been considered directly 
in the calculation. Therefore, how the energy transmission between the longitudinal direction and 
the transverse direction or between the plasma and the air can not be described in detail. In the air, 
they are challenges for our future work. 
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